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ABSTRACT 

The fact that most extrasolar planets found to date are orbiting metal-rich stars lends 
credence to the core accretion mechanism of gas giant planet formation over its com- 
petitor, the disc instability mechanism. However, the core accretion mechanism is not 
refined to the point of explaining orbital parameters such as their unexpected semi- 
major axes and eccentricities. We propose a model, which correlates the metallicity 
of the host star with the original semi-major axis of its most massive planet, prior 
to migration, considering that the core accretion scenario governs giant gas planet 
formation. The model predicts that the optimum regions for planetary formation shift 
inward as stellar metallicity decreases, providing an explanation for the observed ab- 
sence of long period planets in metal-poor stars. We compare our predictions with 
the available data on extrasolar planets for stars with masses similar to the mass of 
the Sun. A fitting procedure produces an estimate of what we define as the Zero Age 
Planetary Orbit (ZAPO) curve as a function of the metallicity of the star. The model 
also hints that the lack of planets circling metal-poor stars may be partly caused by 
an enhanced destruction probability during the migration process, since the planets 
lie initially closer to the central stars. 

Key words: planetary systems: protoplanetary discs - planetary systems: formation 
- stars: abundances 



1 INTRODUCTION 

The discovery of the first extrasolar planet around 51 Pegasi 
(Mayor & Queloz 1995) prompted an intensive search effort 
that continues today. Currently a hundred and sixty of them 
are catalogued, though the definition of planet may not be 
suited to some, whose masses place them more conveniently 
in the brown dwarf realm. The ultimate goal, the detec- 
tion of a planet like our own, lies in a foreseeable future, as 
promised by plans of sophisticated space-based instruments 
capable of detecting Earth mass planets. 

Even with a sample of gas giant planets at hand we face 
more than a few challenges regarding the main competing 
theories of gas giant planet formation, namely, core accretion 
and disc instability. None of them predicted the existence of 
gas giants so close to their parent stars, or the pronounced 
eccentricity of some of their orbits. There is, however, a fea- 
ture of that sample that favours the core accretion mecha- 
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nism: most of the planets are circling metal-rich stars, and 
this relation has been dismissed as a possible outcome of 
selection effects, or the result of pollution by the absorp- 
tion of protoplanetary disc material and/or planets and/or 
disc leftovers such as planetesimals and comets (Santos et 
al. 2003; Pinsonneault, DePoy & Coffee 2001). Santos, Is- 
raelian & Mayor (2004) further argue that the frequency of 
planetary companions is remarkably constant (~ 3%) up to 
solar metallicity, rising very rapidly afterwards and reaching 
~ 30% for twice the solar metallicity, a result that suggests 
the existence of threshold metallicity values for the efficiency 
of planet building. 

Even studies that indicate some influence of selection 
effects and pollution (Murray & Chaboyer 2002; Laughlin 
2000; Israelian et al. 2001) concede that the high metallicity 
is at least partly caused by the primordial metal-rich nebula 
that formed the stars. The fact that the photometric sur- 
vey (Gilliland et al. 2000) of ~ 34,000 stars in the globular 
cluster 47 Tucanae ([Fe/H] = —0.7) has found no transiting 
short-period planets is also an important indicator of the 
link between stellar metallicity and gas giant planet forma- 
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tion. The researchers expected to find ~ 17 planets, based 
on the statistics ol radial velocity surveys. However, the re- 
sult of this survey should not be easily taken as typical for 
low metallicity stars because the very nature of the cluster 
may have created a bias. The high star density may have 
imposed extra restriction regarding planet formation and 
stability, since the radiation from the expected number of 
O stars can be sufficient to remove circumstellar discs be- 
fore planets could form (Armitage 2000). Moreover, close 
encounters between stars in the early stages of the globu- 
lar cluster can also disrupt a disc (Bonnell et al. 2001), or 
eject the planets that were allowed to form (e.g., Sigurdsson 
1992), leading to a population of free-floating planets that 
may or may not remain bound to the cluster, depending on 
the velocity of the ejection (e.g., Hurley & Shara 2002) . 
Altogether, these processes might explain the lack of plan- 
ets in 47 Tuc. It is worth noting, on the other hand, that 
another recent photometric survey of 47 Tuc (Weldrake et 
al. 2005), more sensitive to the uncrowded outer regions, 
indicates that system metallicity is the dominant effect. 

Since both the accretion disc and the parent star are 
products of the same gas and dust cloud, a relationship be- 
tween primordial metallicity and the frequency of planet 
finding would be a natural outcome in the core accretion 
scenario (Pollack et al. 1996), which holds that a rocky nu- 
cleus of around 10 Earth masses is formed before the mas- 
sive infall of gas completes the formation of a gas giant. 
The disc instability model, on the other hand, is based on 
a self-gravitating mass of gas that is formed in the rotating 
accretion disc and is fairly insensitive to metallicity (Boss 
2002a). One of the main merits of the disc instability model 
is the short formation time (10 2 to 10 3 years) of protoplan- 
ets, observed in the simulations by Boss (2002b) and Mayer 
et al. (2002), compared with the considerably longer time 
in a core accretion scenario (up to lOMyr), which would 
in principle violate the age constraint for accretion discs, 
based on observation of young stellar objects (Haisch, Lada 
& Lada 2001; Bary, Weintraub & Kastner 2003). However, 
the simulations do not account for the long-term survival of 
the protoplanets, or the origin of the instability. Moreover, 
simulations incorporating stochastic core migration through 
the protoplanetary disc (Rice & Armitage 2003) show that 
a Jupiter mass planet can be formed within lOMyr. 

If the core accretion mechanism is indeed the one that 
regulates the formation of gas giant planets, as well as of 
rocky ones, then the existence of gas giant planets around 
metal-rich stars is probably but one consequence. 



2 OPTIMUM ENVIRONMENT FOR PLANET 
FORMATION 

Let us consider the key radial variables for planet formation 
in the accretion disc, disc temperature and surface density. 
The disc temperature drops continuously as the distance 
to the star increases, following a power law. For example, 
a crude profile can be derived from gravitational potential 
energy conservation and blackbody radiation (Hartmann 
1998), showing that the disc temperature is proportional to 
r~ 3 / 4 , where r represents the radius; a similar relation can 
be written considering only the radiation from the central 
star that is intercepted by the disc. 



Since the gravitational pull of the embryo gas giant 
planet gets more efficient for low gas temperature, we should 
expect gas giants to form at considerable distances from 
their host stars. On the other hand, the roughly spherical 
form of the contracting nebula dictates that the disc sur- 
face density has a decreasing radial profile (Hartmann 1998; 
Pringle 1981), so that at large distances from the star the 
advantage of low temperatures is hampered by the lack of 
available matter. The "optimum formation radius" is there- 
fore a compromise between these two parameters. We should 
then expect that at this radius the probability of planet for- 
mation is maximum, or, considering multiple planet systems, 
that at this point the most massive planet is expected to be 
formed. 

This scenario is consistent with an embryo planet al- 
ready pulling gas and dust, but the initial step of planet for- 
mation, according to the core accretion model, would involve 
not the surface density, but more precisely the dust surface 
density, since a rocky core is formed in the first place. The 
dust surface density also decreases as distance increases, so 
that we have again an "optimum formation radius" . 

But here we hit a possible connection between the stel- 
lar metallicity and the semi-major axis of our embryo planet; 
for we could then expect that for stars with higher metal- 
licity (and considering the stellar and accretion disc masses 
equal), the radial profile of the dust surface density is al- 
tered. 

According to the studies of the radial flow of dust par- 
ticles in accretion discs, by Takeuchi & Lin (2002) , medium 
size particles (10 to 100 (im) will migrate inward rapidly (or- 
bital decay time ~ 10 6 yr for 100 /im particles), distorting the 
radial profile in a way that in the inner part of the disc (up 
to ~ 15-20 AU), which is the relevant distance range for the 
formation of the class of planets we are discussing, it will 
get flatter (see figure 9a of Takeuchi & Lin 2002), becoming 
steeper further out. Larger particles (^ 1mm) will concen- 
trate in an even narrower region (upper bound ~ 1 AU), be- 
yond which the dust surface density falls to a value around 
one order of magnitude below that of 100 fim particles. 

If we also take into account that gas giants are not 
believed to form through core accretion beyond ~ 20 AU 
(Mayer et al. 2002), since the formation of a solid core would 
take too long with the surface densities involved, the pro- 
file in the region of potential gas giant planet formation is 
flatter compared with that of the outer region, and we can 
infer that it will get even flatter in a higher metallicity disc, 
since more dust particles will migrate radially inward. This 
phenomenon will push the "optimum formation radius" out- 
ward, since the competing effect of a lower temperature will 
tend to get the upper hand. 

Thus, considering the physical variables of both the ac- 
cretion disc and the parent star that were taken into account, 
there may be a significant correlation, in a sufficiently large 
sample of planetary systems, between the initial semi-major 
axis of the most massive gas giant planet in a given planetary 
system and the metallicity of its parent star, since the most 
massive one would probably be formed at or near the opti- 
mum radius. Such behaviour would be constrained at great 
distances, where the dust surface density becomes too low, 
and at very low radii, where the high temperature prevents 
the massive accretion of gas from taking place, leaving only 
small, rocky planets like the Earth. The new planets cir- 
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cling very metal-rich stars and very metal-poor ones would 
then tend to be located at two asymptotic semi-major axis 
(SMA) values. These arguments lead naturally to the sug- 
gestion that a growth curve might exist. 

This model can be parametrized mathematically by the 
following reasoning. 

Let P be the probability of giant gas planet formation as 
a function of radius; it can be assumed as being proportional 
to the dust surface density, as (Lineweaver 2001; Wetherill 
1996), since the planet formation process requires the a pri- 
ori agglomeration of a rocky planetary embryo, whose rate of 
growth is directly proportional to the surface mass density of 
solids (Pollack et al. 1996; De Pater & Lissauer 2001), which 
in turn is directly related to the dust surface density at the 
time of the formation of solid planetesimals. Although the 
results by Pollack et al. 1996 indicate such proportionality 
only during the so called Phase I, which is short compared 
with Phase II, we understand that P depends first of all 
on the formation of a rocky core. Moreover, the time spent 
in Phase I can be substantially increased depending on the 
initial conditions, as shown by Pollack et al. in the case of 
the formation of Uranus. P can also be assumed as being 
inversely proportional to the midplane disc temperature, T, 
since the gas accretion rate that ultimately forms the gas 
giant is directly proportional to the gas density in the im- 
mediate vicinity of the protoplanet (Pollack et al. 1996). 
The gas density, in turn, can be considered proportional to 
1/T, as dictated by the ideal gas law, in a nebula with gas 
in pressure equilibrium, at least until the beginning of the 
runaway gas accretion. If the gas accretion rate is too small, 
the formation time exceeds the lifespan of the gas in the 
protoplanetary disc and no gas giant will be formed. In the 
simulation conducted by Pollack et al. (1996), the tempera- 
ture also influences the gas accretion rate in a more subtle 
way, as explained in the beginning of this section: the effec- 
tive radius of the planet is defined by the boundary beyond 
which the thermal energy of the gas exceeds the gravita- 
tional binding energy, and the gas accretion rate is directly 
proportional to the effective radius, so a planet travelling 
through a cooler gas will tend to grow faster. Here we note 
that the statement by Pollack et al. (1996) that their re- 
sults are insensitive to T and gas density is not valid in the 
context of our approach, since they use pairs of T and gas 
density for the same protoplanetary nebula (Mizuno 1980, 
Bodenheimer & Pollack 1986), in which both variables de- 
crease with radius, whereas we consider that, for a given 
radius, a cooler nebula will correspond to an increase in gas 
density. 

Thus, we can write 



P oc (a s /T). 



(1) 



Considering that the radial profile of as (De Pater & 
Lissauer 2001) is similar to that of surface density as adopted 
in simulations by Takeuchi & Lin (2002) and Boss (2002b), 
then 



as oc r 



(2) 



where a varies between approximately 0.5 and 1.5. 

The disc temperature profile is also assumed to follow 
a power law (Mayer et al. 2002) 

-0 



where in Boss (2002b) (3 ~ 2 and t is a constant (at least 
for a given stellar mass). However, in the more recent sim- 
ulation by Nomura & Millar (2005) j3 ~ 1, which is in good 
agreement with the earlier simulations of DAlessio et al. 
(2001). In Eq. (2) and r is given in arbitrary units. Note 
that even if we consider that the dependence of P on the 
temperature follows a power law, this can be incorporated 
into the parameter j3. 

In order to take into account our hypothesis that the 
probability is also a function of the metallicity, that is 
P — P(r, Z), where Z = [Fe/H], we assume that a = a(Z) 
and /3 = (3{Z), in such a way that the values mentioned 
above (for a and /3) belong to the high metallicity range of 
the functions we are searching for. An increase in Z should 
cause an increase in the dust surface density, lowering the 
value of a, and also an increase in the disc opacity, which 
lowers the temperature gradient. Asymptotic values for a 
and (3 should be assumed when Z is large, since the bound- 
ary conditions ruling the heat exchange process (central star 
photosphere temperature, disc mass of finite extent, and in- 
terstellar matter temperature) always require decreasing ra- 
dial profiles for T and as- On the other hand, large negative 
values of Z must result in large values for a and /3. 

The optimum value of (1) is calculated by setting 

dP dP 
dP = dr + d Z = 0, 
or aZ 

so, dP/dr = and dP/dZ = 0. Using Eq. © and from 
dP/dr = we obtain the optimum formation radius 



1 opt 



(4) 



As expected, a flatter surface density profile will push the 
optimum radius outward. Using dP/dZ — the same pro- 
cedure gives 



r _ 

' opt — 



/d/3 - da\ 
I (da j 



(5) 



so, using Eq. an d @i we arrive at f3 da = a d/3, which 
is valid for a constant ratio between a and /3, that is, j3/a = 
Pi/on, where (3i and oti are, respectively, any inferior limit 
of the sum over some range of the functions a(Z) and (3{Z). 
It is worth noting that even if we consider only dP/dr — 
(and not necessarily dP/dZ = 0) and once it is supposed 
that a successful planetary formation process is the one with 
simultaneous optimum or near-optimum values of a and (3, 
that is 

(3 -a (/3 + d/3) - (a + da) 



t a 



t (a + da) 



(6) 



T oc r~ p +t, 



(3) 



we still obtain j3 da ~ a d/3. 

The relation between a and /3 can be generalized as 

da d/3 

— = -o =-cf{Z) AZ, 
a p 

where c is a positive constant and / is some function of the 
primordial metallicity (which is mirrored by the metallicity 
of the star prior to potential pollution effects). We have put 
the negative sign because of our basic assumption that a (or 
(3) increases as Z decreases. A general form for the so far un- 
known / in Eq. JSJ which takes into account the constraints 
mentioned before and meets the mathematical requirements 
is 
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/(Z) = (l+e«*>) 
where 

C(Z) = c(Z-Z ) 



(7) 



(8) 



The substitution of Z by (Z — Zo) accounts for the fact 
that the value ( = does not correspond necessarily to the 
solar metallicity. Note that / varies from 1 to as Z varies 
respectively from — oo to +oo. Then Eq. 10 states that for 
large values of Z, a and /3 have asymptotic constant values 
and for low values of Z both parameters and their rate of 
change, da and d/3, go to +oo. 

Integrating Eq. JSJ using J7J and @ we have 



dma' = / din/?' 

h 



(l + e^-'dC' 



which gives the functions we have been looking for, respec- 
tively 



a(Z) =a a (l+e-«< z >) 
P(Z)=(3 a (l + e~^) 



(9) 



were a a and /3 a are the asymptotic values at high metallicity 
and £ is given by Eq. JSJ. 

Substituting Eq. © in J1J we have 



' opt 



(Z) = I — I (in arbitrary units). 

V t a a / 



1/(3(2) 



(10) 



Thus, the optimum radius (which can be replaced by the 
SMA of unmigrated planetary orbits) is then linked to the 
metallicity Z. The population of young planets not yet in- 
fluenced by migration will tend to follow a curve, dubbed 
ZAPO, for zero age planetary orbit, in a metallicity versus 
SMA diagram. 

From Eq. 1101 . the minimum and maximum values of 
r opt are, respectively 



lim r opt (Z) 



lim r opt (Z) = 



1, since (3 — > oo 



t tta 



(11) 



To rescale 1 arbitrary unit (in the above equation) to 
AU we take 



1 arbitrary unit = 7 1 AU = 7 1 214.94 R G 



(12) 



so that 7 fixes the most probable minimum distance r opt 
where a massive planet can be formed at a very low metal- 
licitiy. Using Eq. 1101 and 1121 . Eq. @ can be written as 



Tot I - 

.7 



-f3(Z) 



+ 



Pa - «« ( r. 



7 



(13) 



where (r/7) is given in AU. One remark on self-consistency 
is warranted here. The second term of Eq. 11311 is equivalent 
to t in Eq. 0, supposed constant. This second term is actu- 
ally very slightly metallicity dependent, over the metallicity 
interval in which P(r) is non-negligible. As will be shown 
below, in the discussion of the fitting process, P(r) van- 
ishes rapidly for [Fe/H] <— 1, and therefore this second term 
can be considered as constant over the metallicity interval 
of the data. As we describe on the next section, although 
the introduction of the scale factor does not significantly al- 
ter the possible solutions of Eq. 1101 . the slight dependence 



on Z in Eq. 1131 results in temperature profiles for Z down 
to [Fe/H] = — 1 which are more physically meaningful than 
that given by Eq. . 

In order to draw to a crude estimate of a migration pro- 
cess which would displace the entire ZAPO curve to a lower 
value and form different populations of migrated planets, we 
introduce in Eq. 1101 the fraction n, so that 



r(Z,n) 



(14) 



Here we are considering the so-called type II migration 
mechanism, which dictates that the planet (s) formed is (are) 
large enough to open a gap in the protoplanetary disc, mov- 
ing radially in lockstep with the gaseous disc. This migration 
process would eventually stop (Matsuyama et al. 2003), and 
leave the planet(s) in a smaller, but safe, orbit. 

This simple mathematical development, based on op- 
timal conditions, is atemporal by nature, and can not be 
derived from studies of the evolution of planetesimal surface 
density profiles (e.g. Youdin & Shu 2002), because they do 
not provide answers as to when and where the rocky cores 
form. As for the physical credibility of Eq. @, one would 
think that when comparing the surface density a(Z) with 
the temperature profile /3(Z), the latter would tend to have 
a much weaker dependence on the migration of dust parti- 
cles than the surface density. In fact, that a(Z) and /3(Z) 
show the same functional relation should be seen only as 
a rough approximation of physical reality. For example, for 
very low values of Z, (3(Z) would become very large, and 
the disc would become radially isothermal, which does not 
correspond to reality. However, even if we do not consider 
that (3(Z) is optimized in Z, and choose an arbitrary smooth 
function that varies with Z between two finite asymptotic 
values, and substitute it in Eq. (JIJ, we still obtain a curve 
very similar to the one obtained by Eq. 1101 . only that new 
adjustment parameters for the model will appear. Therefore, 
for the purposes of this work, the rather simplistic form of 
/3(Z) and a(Z) is robust enough so that our main conclu- 
sions are not affected. 



3 RESULTS AND DISCUSSION 

We now proceed to fit the model to the observed distribution 
of planetary SMA as a function of stellar metallicity. Since 
we are dealing with the optimum semi-major axis for planet 
formation, we expect that the most massive planet (in each 
planetary system) will tend to be formed at that locus, and 
the Doppler shift technique can discover the most massive 
ones orbiting a given star. There is of course the possibil- 
ity that massive planets stay hidden in long orbits, beyond 
the sensitivity of the technique, or because the time line 
available is not currently long enough to enable detection. 
This is unlikely, since planets with periods reaching ~ 11 yr 
have been uncovered. Their apparent absence in stars more 
metal-poor than [Fe/H] /J, — 0.4 (Santos et al. 2003), however, 
is noteworthy. Against the possible explanation of this ab- 
sence as arising from the larger radial velocity uncertainties 
in metal-poor stars due to their shallower spectral lines, San- 
tos et al. (2003) show that the mean uncertainty in the ra- 
dial velocity determinations of the CORALIE spectrograph 
is only 1-2 m/s higher for metal-poor stars, as compared 
with metal-rich ones. 
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Figure 1. Stellar Metallicity versus Semi-Major Axis (SMA) for 
cxtrasolar planets orbiting stars with masses inside the range of 
1 ± 0.2 solar mass. In the case of multiple systems, only the most 
massive planet was chosen. 
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Figure 2. The same as Fig. 1, taking the range of 1 ± 0.1. 



A possible source of selection effect is the fact that 
very short period orbits are not occupied by massive plan- 
ets (Zucker & Mazeh 2002, Santos et al. 2003), which may 
indicate that the migration process is more efficient for less 
massive planets. However, this phenomenon does not affect 
our main conclusions, as will be seen below. 

In our analysis we consider, for each star that harbours 
more than one planet, only the most massive one and assume 
that there is no significant probability of finding still more 
massive planets. On the other hand, we disregard the "plan- 
ets" with minimum mass above 10 Jupiter masses (Mj), 
since this is likely to be the upper limit for the mass of a 
planet (Jorissen, Mayor & Udry 2001). Besides, if we con- 
sider the critical mass to be 13 Mj, and an aleatory distri- 
bution of orbit inclinations, the mean value of the true mass 
is around 1.3 times the minimum mass, and we have again 
a minimum mass of 10 Mj as a constraint. 



Another point to be considered is the mass range of 
parent stars, from 0.32 to 1.60 solar masses. Our model does 
not depend on the stellar mass, although it influences con- 
siderably the temperature profile of the protoplanetary disc, 
not to mention its mass surface density. Therefore, we sep- 
arate the planet harbouring stars according to mass ranges, 
in order to filter out this influence. Since the search for ex- 
trasolar planets is directed mainly at stars similar to the 
Sun, we limit the analysis to stars whose masses fall within 
two ranges of the solar mass, ±20% (72 stars, Fig. 1) and 
±10% (46 stars, Fig. 2). The former accommodates reason- 
ably well the mass determination uncertainties, which are of 
the order of ~ 0.05 solar masses. As we show in the following 
paragraphs, our conclusions are not affected by the choice 
of this mass interval. 

Data on orbital parameters was extracted from the Ex- 
trasolar Planets Encyclopedia, maintained by Jean Schnei- 
der (Schneider 2005); the metallicity values were taken from 
the literature and from the Extrasolar Planets Encyclopedia. 
Table I contains all data used in this work and the references 
for the data used. 

The stellar age range of the sample thus assembled is 
a significant issue when one considers the timescale of the 
migration mechanisms, which act before the proto-planetary 
disc fully dissipates. We have plotted the stars of Table I in 
the theoretical HR diagrams of Schaerer et al. and references 
therein (1993) and have thus inferred stellar ages. These 
determinations suffer from heterogeneities in the literature 
database of metallicity and effective temperature (for the 
latter, the sources of metallicity were used) , besides the un- 
certainties in the models themselves. Given the fact that the 
ages thus determined have probable uncertainties of a few 
billion years (Gyr), we have divided the sample of Table I 
in stars termed "young" , those with ages inferior to 1-2 Gyr, 
those "middle-aged" like the Sun, with ages within 2 and 6 
Gyr, and the "old" stars, older than 6 Gyr. No trend emerged 
between stellar ages and the planetary orbit parameters. Ac- 
tually, examples of very young stars in our sample are com- 
patible with ages of 0.3 Gyr or less (Porto de Mello & da 
Silva 1997, for HD147513; Kiirster et al. 2000, for HD17051; 
see also Mayor et al. 2004), and it is therefore reasonable to 
assume that migration processes stopped long ago, produc- 
ing multiple populations. However, it is also reasonable to 
consider, given the absence of significant observational con- 
straints for the time being, that the migration mechanism is 
the same for all stars. Although there is not yet consensus on 
how the migration process stops, saving the planets (at least 
the ones that have been observed) from being swallowed by 
their parent stars, a correlation between SMA and metallic- 
ity could still be seen after the migration stops, for the ones 
that suffered a slight migration process. Thus, our model 
seeks to fit the observed distribution of planetary SMA as 
a function of stellar metallicity for the (almost) unmigrated 
planets. For the planets that suffered an intense migration 
process, it would be more difficult to distinguish between 
the different populations. 

Fig. 1 shows the metallicity versus SMA plot for the 
most massive planets orbiting the 72 stars with 1 ± 0.2 solar 
masses. Jupiter is also included (at 5.2 AU), and is assumed 
to be an example of a planet that suffered little or no mi- 
gration (Franklin & Soper 2003). There is a growth curve, 
in light of our model, representing the planets that suffered 
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Table 1. Stellar and planetary data 



Star 


Star 


[Fe/H] 


Ref. 


Planet 


Planet 


Star 


Star 


[Fe/H] 


Ref. 


Planet 


Planet 




Mass 






SMA 


Min. Mass 






Mass 






SMA 


Min. Mass 




M Q 






(au) 


Mj 






M 






(au) 


Mj 


Gl 876 


0.32 


+0.00 


4 


0.21 


1.98 


HD 


73526 


1 


.05 


+0.27 


1 


0.66 


3.0 


HD 114386 


0.54 


-0.08 


1 


1.62 


0.99 


HD 


117618 


1 


.05 


+0.04 
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Data on planet orbits and masses was obtained at http://www.obspm.fr/encycl/catalog.html 

data on stellar metallicity and mass was obtained at the following references: 1= Santos et al. 2004; 2 = Santos et al. 2003; 
3 = Gonzalez et al. 2001; 4 = http://www.obspm.fr/encycl/catalog.html 5 = Bouchy et al. 2004. 
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mild or no migration, and distinct populations of planets 
that have suffered more severe migration. Even if we con- 
strain even more the mass range, to 1 ±0.1 solar masses, the 
growth curve is unaltered, as shown in Fig. 2. 

The fitting process of Eq. 110|l . using Eq. 112|l . produces 
a ZAPO curve (n—1 in Eq. i!4l ) as shown in Fig. 3 where 
the SMA of Jupiter and of the planets orbiting HD 75732 
(5.9 AU) have been considered as constraints (almost unmi- 
grated ones). We should take some value of a a between 0.5 
and 1.5 and /3 a between 1 and 2 (cf. Eq. @ and J3J); the 
constant c in Eq. (|SJ is inversely proportional to some range 
of values for Z where the ZAPO curve grows between the 
asymptotic limits established by Eq. HI IB : Zq is the value for 
the point of inflection of the growing curve. 

The general behaviour of the ZAPO and temperature 
profiles (Eq. tl'Sl ) according to the parameters is the follow- 
ing. Assuming given values of c, /3 and a, and (suitably) 
changing the values of 7, Zo and t, we can fit the data with 
almost the same curve at the high metallicity range while 
the minimum r opt (Z) varies with 7 (at the low metallic- 
ity range). However, low values of 7 are preferred in or- 
der to produce low values of temperature for low values of 
Z and high radius (and thus low values of /3 are also pre- 
ferred). Analogous simulations, now fixing 7 and changing 
c (so changing only the steepness of the curve), show that 
low values of c are preferred to obtain the above behaviour 
of the temperature. 

In order to restrict the choice of the parameters, we have 
done simultaneous fit of the midplane temperature profile 
of Nomura & Millar (2005, their figure 2) (also DAlessio et 
al. 2001, their figure 5a), assuming [Fe/H]= 0, and of the 
ZAPO (n=l), having the two planets cited above as con- 
straints, and also using the requirement that for low values 
of the metallicity the midplane temperature profile should 
be physically reasonable. For a rather good simultaneous fit 
we have a a = 0.8, 0a = 1.0, c = 4, t = 0.028, Z = -0.63 
and 7 = 1.433 (which means that the minimum r op %(Z) 
equals 150 R0, or ~0.7AU). Note that this distance is in 
agreement with the location of the puffed-up inner rim below 
which dust evaporation destroys the dust particles (Isella & 
Natta 2005, Gorti & Hollenbach 2004). As mentioned before, 
for this analysis, it is not necessary to take into account very 
low values of [Fe/H]. We found that, for [Fe/H] < - 1, the 
distribution P(r) becomes very narrow and with low max- 
imum values around values of r below 1AU, such that the 
total probability becomes negligible. 

It is worth noting that, in this formalism, t mathemat- 
ically defines the midplane temperature profile (Eq. ©) in 
a given range of the disc radius, so it does not consider the 
temperature behaviour at the end of the disc. We believe 
that t, a small quantity, may also play the role of the dy- 
namical parameters of the stellar disc, and thus account for 
ZAPO curves for different stellar masses. 

Distinct populations of migrated planets are obtained 
by manipulation of the relation r(Z,n) — nr op t(Z). This 
procedure gives the ZAPO (n—1) and some of the type II 
migration curves (n < 1) on Fig. 3, for the most massive plan- 
ets orbiting stars with masses around one solar mass. We can 
see that down to n = 0.01 there is no clear stratification of 
planet masses, as advanced in the beginning of this section. 

We must also bear in mind that a mild pollution effect 
may have played a part in distorting the curves, particularly 
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Figure 3. ZAPO curve and its evolution due to migration, 
r(Z,n) = nr opt (Z), obtained from Eq. JSJ to {TJ3, and QD< The 

planets were displayed according to three mass ranges, where Mj 
stands for one Jupiter mass. 

for the population of planets with tighter orbits, since their 
central stars must have engulfed more material (gas, dust 
and planets) in order to pull the (remaining) planets that 
close. Evidence that this possible offset in metallicity is lower 
than 0.1 dex in [Fe/H] is provided by Santos et al. (2003) 
and Pinsonneault et al. (2001), who show the absence of any 
trend of possible metallicity enhancements with the mass 
of stellar convective envelope. An effect of this magnitude 
would only displace the stars slightly to the right, without 
altering the global shape of the data points and the curves. 



4 CONCLUSIONS 

The sample of extrasolar planets around stars with masses 
similar to the Sun produces a distinct population of mi- 
grated planets. However, the ZAPO curve, which is all that 
we hope to model so far, is necessarily a tentative one, since 
we have a scarce number of planets with larger SMA. This is 
probably because the Doppler shift technique requires more 
observational time in order to pinpoint larger orbits, and we 
expect that in the near future the ZAPO curve will be more 
easily delineated. It is also supposed, in light of the core 
accretion mechanism, that the larger orbits will not exceed 
much the current maximum, since the steep decline in dust 
surface density would require too much time for planets to 
form in the outer regions of the protoplanetary disc. The 
disc instability simulations, on the other hand (Mayer et al. 
2002; Boss 2002b), indicate that large planets could form at 
even larger distances (> 10 AU). 

The scenario of our model predicts that the metal- 
poorer stars have fewer giant planets not only because the 
dust surface density profile is affected, but also because the 
ones formed lie initially closer to the central stars, and have 
therefore been more frequently absorbed by the central stars 
as a consequence of the migration process. Given the statis- 
tical significance of the lack of long period planets in metal- 
poor stars, this is a prediction of our model that may be 
borne out by future observations, as the sample of detected 
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cxtrasolar planets steadily increases. Also, we note that an 
enhanced probability of planet engulfment for the metal- 
poor stars would reflect in an enhanced efficiency of effec- 
tively observed planetary companions in metal-rich stars. 
The distribution of planet hosts would qualitatively be a 
steepened version of the ZAPO curves, since for metal rich 
stars the migration process would tend to leave behind mul- 
tiple populations of survivors, but for metal-poor stars the 
close-in optimum formation radius would make for a more 
efficient engulfment of planets, leaving behind an observa- 
tional void. Qualitatively, this is in very good agreement 
with the results of Santos et al. (2004; see their figure 7) 
and also explains recent statistical studies by Sozetti (2004) 
indicating that close-in planets are more likely to be found 
around metal-rich stars. Therefore, the arguments used to 
put forward the ZAPO model can explain the observed fre- 
quency of planetary companions with metallicity, but the 
argument cannot be turned around to obtain the former 
from the latter. Moreover, the model provides a physical ex- 
planation for the lack of long period planets in metal-poor 
stars. 

We have to emphasise that the model presented by 
this work is very simple, and that undoubtedly other fac- 
tors should enter the picture in the future, as more obser- 
vations become available. For instance, the quantity P, that 
we called "probability of giant planet formation" is likely to 
be related to the total density of solids, not only dust. It 
is a known fact that the density of solids in a disc meets a 
maximum where water vapour condenses into water ice. At 
this distance - the snowline -, the most abundant volatile 
becomes solid and its contribution dominates over that of 
rock. For a disc of solar abundance, the presence of ice ex- 
terior to the snowline increases solid surface density by a 
factor of four. Stevenson & Lunine (1988) propose that dif- 
fusive redistribution of water vapour through the snowline 
further enhances the density of solids by 1 to 2 orders of 
magnitude, naturally leading to the emergence of a "best 
radius" for giant planet formation. In this sense, a/T can 
be regarded as an approximation to the exact formulation 
Cdust + o"iccs , with 1/T reproducing the process of condensa- 
tion of water vapour into ice as the temperature decreases. 
It is interesting to note that although the discontinuous be- 
haviour of the condensation does not emerge from this linear 
approximation, the location of the peak of our probability 
profiles (for solar and higher metallicity) reasonably agrees 
with the location where the density of solids peaks when ice 
condensation is taken into account (Kornet et al. 2004). 

We acknowledge the fact that the parameter n, which 
embodies our ignorance on the migration process, will prob- 
ably turn out to be quite complex. The very fitting process 
of the curves shown in Fig. 3 suggests that the so called 
type II migration process is not the only one that influenced 
the final orbits. However, even with the approximations as- 
sumed in this work, we believe that the model, coupled with 
a metallicity versus SMA plot, is a simple and useful tool 
for probing the evolution of planetary populations. 

The model predicts the optimum region for planet for- 
mation, but offers no clue as to the mechanism that triggers 
the process. A possible answer lies in the theory that the gen- 
eral radial gas density profile is probably subject to many 
local density enhancements, which attract solids due to pres- 
sure gradients and trigger planet formation (Haghighipour 



& Boss 2003). Therefore, the density enhancement located 
near or at the optimum point would naturally evolve more 
rapidly. 

Finally, it is also worth mentioning that the same model 
presented in this work can be tested against the formation 
of rocky planets, as soon as sizable samples of these become 
available, only that this time the temperature factor shifts to 
the temperature that allows the existence of dust particles. 
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